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A unified model of electron penetration into linear plasma wakefields is formulated and 
studied. The optimum angle for side injection of electrons is found. At smaller angles, 
all electrons are reflected radially. At larger angles, electrons enter the wakefield with 
superfluous transverse momentum that is unfavorable for trapping. Separation of incident 
electrons into penetrated and reflected fractions occur in the outer region of the wakefield 
at some "reflection" radius that depends on the electron energy. 

PACS codes: 



1. Introduction 

Active studies of p lasma wakefi eld acceleration initiated by Taiima fc Dawson ( 19791) 
for laser drivers and Chen et MI ([1985) for electron ones have reached a point where 
quality of accelerated beams comes to the fore. Various injection techniques are now 
under investigation in search of relia ble ways of p r oduci ng low-emittance low-energy 
spread electron beams, as reviewed bv lEsarev et al. ( 20091 ). One promising technique is 
side injection of a low energy electron bunch to the w akefield at some small angle. This 
method was first proposed bv iKalmvkov et oZI ( 2006[ ) to increase trapping efficiency by 
taking advantag e of a specific wave nonhnear ity. Later, injection at an angle was theoret- 
ically studied bv lLuttikhof et al\ (|2007ll2009h as a cure for ponderomotive scattering and 
deleterious effects on the vacuum-plasma transitio n. Recently, side inj ection has found 
application in proton driven wakefield acceleration (|Pukhov et a/.ll201l[ ) as the means to 
save electrons from scattering at early stages of driver evolution in the plasma. 

The theory of side injection is still a long way from co mpletion. It is ess e ntiall y two- 
dimensional. The one-dimensional Hamiltonian approach of Esarev fc Pillo^ ( 1995[) helps 
to some extent to identify rough scalings, but quantitative results necessarily rely on 
simulations. In this paper we make a step toward creation of a two-dimensional trapping 
theory and formulate necessary conditions for side injected electrons to penetrate into 
the wakefield. We assume electrons are injected at a small angle a to the direction of 
driver propagation and show that there is an optimum value of this angle. 

We first formulate the necessary mathematical framework for particle beam-driven 
wakes (section [2]), as there is a universal field asymptotic in this case. In section |3] we 
analyze the equation obtained and find unified laws of electron penetration into the 
wakefield. In section |4] we generalize the formalism to laser-driven wakes. 
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Figure 1. The geometry of the problem in laboratory (a) and co-moving (b) coordinates. 



2. Particle beam-driven wakes 

We use cylindrical coordinates (r, Lp, z) with the drive beam propagating along the 
z axis. The wakefield is assumed axisymmetric and unchanged in time in the driver 
reference frame. This assumption is justified if electron trapping occurs at a much shorter 
time scale than the time of driver evolution, which is typically the case. The phase velocity 
of the wakefield is 14, « c, the relativistic factor is F^, = (1 — / c^)~^/'^ , where c is the 
speed of light. 

The geometry of the problem is shown in figure [1] Electrons initially form a stream of 
the velocity V directed to axis. We use the co-moving coordinate ^ = z — V^t and do not 
Lorentz transform the fields into the moving frame. In new coordinates, electrons move 
backward and radially with velocity components 



V cos a — Vw 



-V sma 



-ac. 



(2.1) 



As far as the side injection is usually aimed at further acceleration of electrons, we assume 
Vuj > V and < 0. 

Electron energies discussed in the context of side injecti on mostly fall in the range 
3-15 MeV (jLuttikhof et aLll2007l [2009t IPukhov et aLll201l[ ). and usually the condition 
a < F~^ is fulfilled, where F = (1 — V"^ /c^)~^^'^ . As we show below, for the optimum 
angle the stronger condition aF ^ 1 must be fulfilled, and we use this to neglect the 
angular dependence of the longitudinal velocity: 
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2F2 / 
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2F2 J 



V-V^ 



(2.2) 



The force exerted on an axially moving relativistic electron is the gradient of the 
wakefield potential energy $. Side injected electrons either fall into the potential well 
of the wakefield, or are refiected radially by a potential hump. The initial energy Wr of 
electron transverse motion is determined by the radial momentum p,, w aVmc, where m 
is the electron mass, 

W. = ^^ —mc^. (2.3) 

For small a, this energy contains the product of two small parameters (a and aF) and 
is much smaller than height of any potential hump in wakefield structures of interest. 
Consequently, the choice of whether electron enters the wakefield or is reflected radially is 
made at large radii. Once an electron is trapped by the potential well, its radial velocity 
quickly increases, and the electron approaches the axis with little change in ^-coordinate. 
Thus we are interested in potential behavior at large radii where important characteristics 
of trapping are determined. 

For the particle driver of charge density p{r,£^) — Pbf{r)g{£,) and the linearly respond- 
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ing plasma, the wakefield potential energy is (|Cheijll987h 

Airpbe 



'^'(^'0 = (2.4) 



1.2 



where 

R{r)^kll f{r')Io{kpr')Ko{kpr)r'dr' + klj f{r')Io{kpr)Ko{kpr')r'dr', (2.5) 
ZiO = kp / sin [kpi^ - 0] , (2.6) 



C V TO 

uq is the plasma density, e > is the elementary charge, Iq and Kq are zeroth order 
modified Bessel functions. At large radii and behind the driver, formulae (|2.4I) - (|2.6I) take 
the universal form 

<I>(r, = $0 cos{kp^ + MKoikpr) (2.7) 

with the amplitude $o a-nd phase 0o determined by the individual driver shape. Setting 
(/>o = by the choice of origin and using the asymptotic form of Bessel function Kq, we 
simplify the potential energy to 

j —k T 

$(r, « <f cos(fcpa^f (2.8) 

Separation of incident electrons into trapped and reflected fractions occurs at some radius 
To ^ k~^ which we define later. Denote $i — $(ro,0). Then, in the vicinity of tq, we 
can write 

$(r,,f) « $icos(fcp0e~''''^''~''°^ (2.9) 
Equations of electron motion in the potential (|2.9I) are 



^ = -^ = kp^i cos(fc,0 e-'''^^-'-o\ (2.10) 



dr pr d£_ 

Here we neglect the change of total electron energy since it happens on much longer 
timescales. Equations (|2.10p - (|2.1ip can be combined into one: 

^ -J^^os[k,iv,t + io)]e-'^^^-^''\ (2.12) 
Introducing dimensionless variables 



i=-kpv^t, |o = fcp6, ■f=|o-i, x = kp{r - tq) - In ( ^^^^^ ] , (2.13) 



we rewrite equation ()2.12p in a universal form: 



^=cos(i-eo)e-^ (2.14) 
This equation must be solved with initial conditions corresponding to electron arrival 
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Figure 2. Family of electron trajectories for v — —0.7 (a) and v — —1 (b). Lower graphs show 
the location of potential wells and humps. 
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Figure 3. Area of possible electron trapping on the plane (^t, v) (a) and total penetration 
efficiency rjt as the function of the initial velocity v (b). The lower graph shows the location of 
potential wells and humps. 



from large x with the negative dimensionless velocity 

dx Vr ac 



V = 



dt V - 



(2.15) 



3. Penetration of electrons into the wakefield 

Equation (|2.14|) can be easily solved numerically for any values of initial velocity v 
and phase ^o- For low initial velocities, all electrons are reflected by the outer region 
of the wakefield [figure [2ja)]. For larger w, there are two types of electron trajectories 
[figure [2l^b)]. Depending on the phase, electrons either penetrate the wakefield, or are 
reflected. The electrons entered into the wakefield are quickly accelerated radially and 
stick to some phase ^t. This is exactly the phase of the wakefield into which the injection 
of electrons occurs. 

General picture of electron penetration is shown in figure |3l[a) . Electrons cannot enter 
the wakefield at phases and initial velocities marked by the white color. The darker the 
color the greater the flux density Ft of entered electrons, i.e. the number of electrons 
trapped in the unit interval of divided by the number of incoming electrons passing 
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Figure 4. Flux density Ft of trapped electrons as the function of final phase for different 
values of the initial velocity v. Thin line is the flux density of incoming electrons. The lower 
graph shows the location of potential wells and humps. 



through a unit interval of ^ at large radii. The total penetration efficiency rjt (into all 
phases) is shown in figure [Sjb). The penetration threshold is observed at w w —0.81. 
For > 1, not only more than half the electrons penetrate into the wakefield, but also 
the electron flux density increases several times (figure S]). In other words, injection of 
electrons preferably occurs into a certain phase of the wakefield, and this phase depends 
on the dimensionless radial electron velocity v. 

Penetration of an electron into the wakefield does not mean this electron is trapped 
by the wakefield. The initial radial velocity of the electron can be too high, or the energy 
gained by the electron in near-axis regions can be insufficient for trapping. In these cases, 
the electron crosses the axis and leaves the wakefield on the opposite side. The lower the 
initial radial velocity of the electron, the easier the electron to be trapped. Consequently, 
there is an optimum angle of electron injection corresponding to \v\ ~ 1. It is just above 
the penetration threshold: 

aopt ~ ■ (o-lj 

c 

If the relativistic factor of the wave is much greater that that of the electron, F^ ^ F, 
then we obtain the engineering formula 



1 



"opt 



2F2 



(3.2) 



For Tm ^ F, the optimum angle is even smaller. Estimate p.2p justifies the condition 
aF < 1 and formula (|2?2|) . 

As we see from figure [2l separation of incident electrons occurs at i « 0. Let us find 
the corresponding dimensional radius vq. Substituting i = and r = tq into definition 
of X (|2l3l) . we find 

$(ro,0) =Fmw|. (3.3) 

Using approximation (|2.8p for $i, assuming F^, ^ F, and limiting ourselves to the 
logarithmic accuracy, we obtain a quick estimate for the reflection radius: 



fcpMn 



(3.4) 
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4. Laser beam-driven wakes 

For laser drivers, there is no universal asymptotic form of the wakefield potential at 
large radi i. The wakefield po tential energy is determined by the vector potential of the 
driver A (jEsarev et a/.ll2009l ): 

Hr,0 = -mc'kpj^ -A^sm[kp{^' -0]de, a=— . (4.1) 

Usually $(r, ^) follows the same radial dependence as A^(r, ^) does. This dependence 
does not necessarily have an exponential asymptotic at large radii. Nevertheless, the 
results of section [3] are still applicable to laser driven fields though with a lower accuracy. 
We can approximate the potential energy in the vicinity of separation radius tq by the 
proper exponent function: 

$(r, « *i cos{kp£,) e-^^'-'^o) . (4.2) 

For example, for a Gaussian driver of the radius ar we have 

cx e-'-'/"', $(ro + Sr,^) = $i cos(fcpO e-^-o^rV-^-^r-V-^ ^ x = 2ro/al (4.3) 

As far as A 7^ fcp, equation (|2.12p . definition of x (|2.13p . and the initial velocity (|2.15p 
need to be modified: 

^ = ^cos[kp{v,t + ^o)]e-'^^-^"\ (4.4) 

--Hr-ro)-ln(^-^^y = | = (4.5) 

In newly defined variables, equation (|2.14p remains unchanged, as are the results of 
section |3] presented in figures [2H11 We need to modify only expression p.ip for aopt and 
equation (13.31) for r^: 

aopt , ''^'°^'^^^''^ = ^^)[-^^) (4-6) 



5. Summary 

Let us summarize the main findings. There is an optimum angle for side injection of 
electrons into the plasma wakefield, as is given by formulae p.ip or (j4.6p . At smaller 
angles, all electrons are reflected back. At larger angles, electrons enter the wakefield 
with superfluous transverse momentum that is unfavorable for trapping. 

Penetration efficiency and the wakefield phase into which electrons enter are deter- 
mined by the field behavior at large radii. In most cases, the wakefield there decays 
exponentially, and the process of electron entrance into the wakefield can be studied in 
a unified way. 

Although formulated for linear wakefields, the results of the paper are also applicable 
to nonlinear wakes if the latter have a factorable field asymptotic of the form (|4.2p . 

This work is supported by RFBR Grants 09-02-00594, 11-02-00563, and by the Russian 
Ministry of Education grants 2.1.1/3983 and 14.740.11.0053. 
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